arXiv: 1509.01491 vl [physics.class-ph] 3 Sep 2015 



DICAUS 


Dr.Paul.Kinsler@physics.org 
Causal Diagrams http;//www.kinsler.org/physics/ 


Causal diagrams for physical models 

Paul Kinslef] 

Blackett Laboratory, Imperial College London, Prince Consort Road, London SW7 2AZ, United Kingdom. 

(Dated: Monday 7* September, 2015) 

I present a scheme of drawing causal diagrams based on physically motivated mathematical models expressed 
in terms of temporal differential equations. They provide a means of better understanding the processes and 
causal relationships contained within such systems. 


I. INTRODUCTION 


This paper presents a systematic way of making diagram¬ 
matic representations of causal physical models. The premise 
follows that due to the presentation of Kinsler IT], with a no¬ 
tion of causality being linked to the existence of temporal 
derivatives in physical models based on differential equations. 
Whilst many other approaches to causality are both useful and 
relevant to many areas, the one here is chosen to be simple and 
direct: it analyses the causal properties of an existing model, 
and acts as a guide for the constmction of models consistent 
with its principles. 

The original motivation for these diagrams was in analogy 
with the fishbone-like Ishikawa diagrams 0, used in busi¬ 
ness, and created by Kaora Ishikawa in order to identify the 
origins of a particular effect. However, unlike in business 
where the analysis, no matter how carefully undertaken, is far 
from being completely specified, in physical diagrams I can 
add features to precisely represent the physical model. Note 
that a simple way to “physicsify” Ishikawa diagrams had been 
implemented by Thomas Wong 0, in a third year undergrad¬ 
uate essay project at Imperial College London which I pro¬ 
posed and then supervised in 2015. 

The diagramming process introduced here is distinct from 
that used in standard block diagrams i). In block dia¬ 
grams, elements have specified and defined inputs and out¬ 
puts; whereas in my scheme each node indicates what modi¬ 
fies the quantity of interest. Further, block diagrams empha¬ 
sizes solutions (e.g. by representing integrators) or frequency 
domain properties; whereas my aim is to emphasize on-going 
dynamical behaviour. In particular, block diagrams often loop 
outputs back to inputs to represent feedback, whereas here 
such loop-like constructions would represent non-causal be¬ 
havior. 

Another type of diagrams used in causal analyses are “di¬ 
rected acyclic graphs” (DAGs), and in addition to more tech¬ 
nical aspects m have been proposed as a methodology to as¬ 
sist draw causal inferences from non-experimental data (see 
e g- 0 ). When applied in the contexts for which they are in¬ 
tended, which is not usually for physical modelling, they are a 
valuable tool for investigating and clarifying causal relation¬ 
ships. One of the promotional points about DAG’s is that they 
are a mathematics-free method that relies on graphical rules 
for their constmction. However, as a theoretical physicist, I 
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would prefer any graphical rules to also be a re-representation 
of a specific mathematical model - and this indeed is the case 
for the diagrams I propose here. Thus my proposed diagrams 
do not only say “this causes an effect on that”, but also contain 
a mathematical specification of it. Nevertheless, one might 
start by constructing a diagram for a causal physical model 
in a DAG-like way, based on expectations of cause and ef¬ 
fect. Then, subsequently, this could be augmented with the 
diagrammatic notation introduced here to specify it more ex¬ 
actly; although note that DAGs whilst can be (and are) explic¬ 
itly associated with a labelling of cause and effect, they do not 
have the same sense of temporal dynamics as do many physi¬ 
cal theories. For example, a causal DAG has no way, without 
an augmentation such as the addition of temporal labelling, of 
diagramming an oscillatory system. 

In section [I^ I present the basic core diagrammatic lan¬ 
guage, consisting of arrows and dots and how they relate to 
simple temporal differential equations. This is followed in 
section[In]by extensions allowing multipliers, parameters, and 
other operators to be added to the diagrams by adding labelled 
boxes, as well as how quantities that affect themselves should 
be shown. In section IV I add an optional arrowhead notation 
to simplify the inclusion of spatial derivatives, and remark on 
how simple second order wave equations can be drawn. After 
that, in section |V] I show and discuss interconnected systems 
such as wave models, and in particular how the macroscopic 
Maxwell’s equations can be diagrammed. Finally, I conclude 
in section rvm 


II. CORE DIAGRAMS 


As defined by Kinsler IT], the simplest possible causal 
model for the behaviour of some quantity R under the influ¬ 
ence of a stimulus 2(f), is the temporal differential equation 

‘^ = d,R = Q{t), (2.1) 

dt 

although we might also straightforwardly generalize this to 
models with n-th order time derivatives as 

^=dpR = Q{t). ( 2 . 2 ) 

Note that even though here Q(t) is some function of all time, 
and might even have been defined in advance in some pre¬ 
determined way, the current value of R is only determined by 
the current or past values of Q. In more complicated models. 


1 









DICAUS 


Dr.Paul.Kinsler@physics.org 
Causal Diagrams http://www.kinsler.org/physics/ 



FIG. 1: Primitive vertices and their notation: simple equations with 
only a first or a second order time derivative. 


it may also be that the stimulus Q is (or is calculated from) 
some other quantity which has its own dynamical behaviour. 

For these simple cases, the core diagrams denote each 
causal relationship by means of a simple vertex with a dot 
or filled circle. We denote any quantities that are evolving in 
time, i.e. that are being affected with horizontal lines going 
from left to right. Other stimuli (causes) that affect it come in 
as diagonal arrows pointing either up-right or down-right. 

In fig. [^we show the most primitive vertex with a known 
stimulus Q driving first-order temporal changes in R. If the 
model has higher order time derivatives, the “dot” vertex is 
augmented by one extra circle per extra time derivative. 

It is important to note that the time derivative or deriva¬ 
tives are only denoted by the dot and its surrounding circles 
(if any), and not by the arrowhead. In fact, from a purely 
diagrammatic perspective we could dispense with the arrow¬ 
head entirely. However it helps reinforce the intent that time 
is advancing as we move from left to right across the page, 
and that some stimuli (written Q in the above) is applied and 
drives changes in a quantity of interest (R above). 


III. SCALINGS AND COMBINATIONS 


Sometimes - or indeed typically - the various “cause” stim¬ 
uli we might use are scaled or changed versions of other quan¬ 
tities or stimuli in the model, so on fig. [^we show how to di¬ 
agram such cases. Note that the box notation used here could 
also be used instead of the varied arrowheads suggested below 
as a shorthand for common spatial derivative terms. 

For a a system influenced by two different stimuli, e.g. S 
and Q, the causal diagram is shown on fig. where we just 
put the mathematical terms in a box where they can be easily 
read and understood. We might say that these boxed quantities 
tell us how the stimuli/causes S,Q are “conditioned” before 
altering the behavior of the quantity R. While these condi¬ 
tioning (boxed) terms are, in terms of our model, completely 
specified, the physical processes they represent are very likely 
to be very complicated. If we were to develop a causal model 
that explained the /(2) term in fig. for example, it may well 
generate a complicated interconnected mesh-like diagram in¬ 
volving a variety of relevant quantities and/or stimuli. 

Note that to avoid ambiguity it is necessary to put the full 
expression in each conditioning box. For example, in fig. 
we must write —in the box(es) conditioning S. If we were 
to put only —/r and assumed that a multiplication by S was 
implied, an alternate conditioning term —fiU^ derived from a 


FIG. 2: Scaling and combination of causes. The left hand diagram 
shows the two causes S and Q with their rescaling or functional de¬ 
pendence indicated in their respective boxes each on different arrows. 
The centre diagram shows the same situation as the left, but with the 
modified causes —jlS and f(Q) added together (shown by the joining 
of the the two lines) before application to R as indicated by the ordi¬ 
nary arrowhead. The right hand diagram again re-represents the left 
hand one, but also introduces the diamond element as a way of intro¬ 
ducing a combined stimulus T. Although not strictly necessary here, 
it is an option which can be useful in more complicated diagrams. 


stimulus U would require us to put the potentially confusing 
—jdU in the box. 






dtR = —yR + F d^R = F — yd,R— cOqR 


FIG. 3: Left: Simple linear decay, with a driving term F and a self¬ 
acting term representing a loss proportional to y. Right: Simple 
Lorentz oscillator model, with a driving term F and a self-acting 
term representing a loss proportional to y, and a natural resonance 
frequency cOq. Where each arrow terminates is important: the loss 
term is a first derivative effect and so its arrow ends on the central 
dot; the driving and resonant frequency terms are second derivative 
and terminate on the circle. 

Many quantities of interest have behaviour that depends on 
their current state - a simple driven, linearly damped system is 
one such. For a loss parameter y and driving term F, the causal 
diagram is shown on the left hand side of fig. The Lorentz 
oscillator model in fig. |^is a little more complicated. It has 
two parts - a first derivative “loss” part with dfR —yR needs 
to be added to a second derivative “oscillation & driving” part 
with d^R —(OqRFF. To do this we turn the first part also 
into a second derivative; i.e. —> d^R oc —ydfR, to give the 
expected (Lorentz) sum, 

d^R = -yd,R - (oiR + F. (3.1) 

The outer circle, denoting the second derivative effects, is here 
drawn much larger than in fig. because of the need to vi¬ 
sually distinguish between the effect of the loss term whose 
net effect is first order (~ dfR) and the oscillation and driving 
terms whose net effects are second order (^ d^R). 

Other systems we might make diagrams for are simple 
waves. In a scalar wave, the temporal response of the wave 
amplitude depends on the spatial changes in profile of the 
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wave. Two common cases are ones with first order time 
derivatives, such as the Schrodinger wave equation, and those 
with second order time derivatives, such as the Helmholtz 
equation. Shorn of any parameters such as wave speeds in 
order to simplify the diagrams, these are shown on fig. 


R R 



d,R = VVR d^R = VVR 


FIG. 4: Simple single-field second order wave equations, where a 
scalar field quantity R acts on itself after conditioning: variation in 
its spatial profile causes changes in the temporal evolution. On the 
left hand side we have diagrammed a diffusion or Schrodinger-like 
equation, with their characteristic single time derivatives, whereas on 
the right hand side we have a Helmholtz-like equation. As usual, the 
spatial gradient operator is defined with V = (dx,dy,d^). 


IV. SPATIAL DERIVATIVES 

As we saw in the latter part of the previous section, and 
indeed in many interesting physical systems, the “cause” af¬ 
fecting the quantity of interest is a spatial derivative of some 
kind. It is therefore useful to utilize a shorthand notation of 
different arrowhead types to denote the three vector calculus 
spatial derivatives of interest, which we denote by using dif¬ 
ferent kinds of arrowhead, as shown on fig. 


R fi R 



4R = Ve 3,R = VQ a,R = VxQ 


FIG. 5: Simple vertices involving spatial derivatives; with an open 
arrowhead for the usual gradient (V), a closed arrowhead for the di¬ 
vergence (V •), and a diamond arrowhead for the curl (V x). Here, a 
boldface symbol denotes a vector quantity. 

For multiple derivatives, we can simple stack the arrow¬ 
heads up in the correct order, which is the order in which the 
operator they represent is applied to Q, as seen on fig. I^and 
fig- 0 Further, we can re-diagram the wave models of fig. 
1^ The result of this process is shown on fig. where sim¬ 
ple parameters, a wave speed c and a length scale A are also 
included. Note the ordering of the box and arrowheads in re¬ 
lation to the equation. As we follow the “cause” line, we meet 
the boxes first (so that e.g. R c^R), and then the derivative 
arrowheads (so that c^R —> V • V {c^R)). 

At this point we have all the tools necessary to draw causal 
diagrams of any physical model, at least in principle. How¬ 
ever, in the next section I will apply the concept to some more 
specific models in order to see how the method works out in 


fi R R 

3,R = yyQ (9,R = VxVxQ 5,R = VVQ 


FIG. 6: Vertices involving double spatial derivatives. Note the or¬ 
dering of the arrowheads, which may at first seem counter-intuitive. 
The arrowheads appear along the line towards the dot in the order 
with which they are applied to Q, thus in the left hand diagram the 
open arrowhead (for V) appears first, then the closed arrowhead (for 
V •) second. The reason for this can be seen in the right hand diagram 
on fig.[^ 


R R 



d^R = y-yQ 52r = vxVxq 4-r = vvq 


FIG. 7: Vertices involving double spatial derivatives driving doubled 
time derivatives. 


practice. One point to note is that is can be tricky, especially 
in models with multiple quantities being affected (“effected”), 
to keep the diagrams well organized enough to remain read¬ 
able. While a valuable process in itself, because it forces one 
to be very clear about cause and effect as implemented in the 
mathematical model, it may be that restricting the use of these 
diagrams to fragments of a larger model is their most practi¬ 
cal use. The challenge results largely from the fact that each 
quantity needs to be represented by a line that represents its 
a temporal history. In contrast, non-dynamical diagramming 
schemes like DAGs or block diagrams, which represent rel¬ 
evant quantities as discrete localized elements such as letters 
or blocks, allow greater freedom when drawing multiple in¬ 
terconnections. 


V. INTERCONNECTED PROCESSES 

Many physical systems have more than a single quantity 
undergoing changes due to their environment; thus, in our di¬ 
agrams, we may need more than one horizontal line. One no¬ 
table type of interconnected system is a wave, where (e.g.) the 
displacement and speed profiles of a stretched string are cou¬ 
pled together in a way that forms oscillations that travel along 
it. Waves are therefore an ideal test case for the diagramming 
of interconnected systems. 

However, let us first consider a simple oscillatory system, 
where the displacement x of some object changes according 
to its velocity v, and its velocity v changes according to some 
(restoring) acceleration which is proportional to its displace¬ 
ment. This is shown on fig. along with an alternate repre¬ 
sentation based instead on displacement x and momentum p, 
with parameters mass M and restoring force constant K. 

In the case of causal wave models, we find that spatial varia¬ 
tion in the wave amplitude drives (causes) changes in the wave 
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R R R 



3,R = VV{AcR) S^R = VV(c^R) 3fR = Vc-VR 


FIG. 8: The wave diagrams for fig. but using the new arrowhead 
notations. Also added are some necessary wave parameters, a length 
scale A and a wave speed c. In the right hand diagram we can see the 
result of moving the parameter in between the spatial derivative 
terms, and by comparing with the centre diagram, also see how the 
arrowhead ordering has been preserved. Note that the open arrow¬ 
head denoting the gradient V might now be considered redundant, 
since it must also be included in the boxed term; however it does aid 
the readability of the diagram. 

d,x = v d,x = plM 

: WT : wr 

d,v = —F^x d,p = —Kx 


FIG. 9: A diagram for a simple oscillating system such as a mass 
on a spring or a pendulum. The left hand diagram shows the sys¬ 
tem based on displacement x and velocity v, where the oscillation 
frequency depends on F\ the right hand diagram shows the system 
based on displacement x and momentum p, where the oscillation fre¬ 
quency depends on M and K. 


- as enclosed by the dotted ovals - is needed to represent the 
model. 


d,P = V -v 



dfy = c^VP 


FIG. 10: A simplified p-Acoustic wave equation, i.e. a sourceless 
scalar-l-vector wave equation, consisting of two first-order parts. We 
leave the wave-speed squared “conditioning” parameter off the 
diagram (but not the equations) for clarity. A unit cell is indicated by 
the dotted oval. 


^(E = £-^(VxH + J) 



profile m. Although second order wave equations, as shown 
in the previous sections, can be diagrammed rather simply 
with a single horizontal line, these hide hide half of the dy¬ 
namics involved in the wave. From some points of view |(7I 
wave equations exhibiting second order derivatives are merely 
contracted forms based on a pair of first-order equations, such 
as in EM or p-Acoustics H. Such models have two evolv¬ 
ing fields, and hence a pair of horizontal right-travelling lines, 
where each fields is driven by the spatial derivative of the 
other. 


Perhaps the simplest example of this is p-Acoustics 0, a 
scalar wave theory usually represented by two scalars and two 
vectors linked by first order differential equations and consti¬ 
tutive relations. To start by keeping things simple, however, 
we can reduce the model to only two quantities, a scalar pop¬ 
ulation P distribution and a velocity field v. On fig. [TO] we 
show this simplified p-Acoustic model with trivial constitutive 
properties (i.e. only a single speed c). In this model the spatial 
divergence in the “flow” vector field v cause the “occupation” 
scalar field P to change in time; whilst spatial gradients in the 
occupation field P cause the flow field v to change in time. 
Likewise on fig. 11 we do the same for a simplified EM, for 
the electric field E and magnetic field H; we see the ordinary 
curl vector equations as influenced by a current source J, with 
only a scalar permittivity e and permeability jj.. 


Eor both figures [T0| and [TT] only a small segment of the to¬ 
tal time history that would otherwise repeat endlessly into the 
future is shown. Strictly, only a single subunit or “unit cell” 


FIG. 11: The diagram for the vacuum Maxwell’s equations, which 
consist of two first-order equations in two the vector fields E and 
B. We leave the constitutive parameters off the diagram (hut not the 
equations) for clarity. A unit cell is indicated hy the dotted oval. 


A more complete causal diagram for Maxwell’s equa¬ 
tions is shown on fig. It mimics the DPBM-hased 

EDTD scheme for their numerical solution, and allowing for 
a Lorentzian response for both dielectric polarization P and 
magnetization M. Although this representation does not need 
to include the fields E and H, since they follow directly from 
the in-diagram known D/eo — P and B//ro — M/fio, it can 
helpful to add them to the diagram to aid readability, as I have 
done. 

In fact, the link with DPBM numerical scheme is not an ac¬ 
cident - both these diagrams and a computer simulation of a 
dynamical system necessarily must be structured in a causal 
way. Computer programs have to calculate things in compu¬ 
tational order, and since this is most likely to be in time order, 
the demands of these diagrams and simulation are the same. 
Note, however, that not all dynamical simulations need to 
chose the time axis to solve along - it is common in optics, and 
sometimes in acoustics, to choose a spatial axis along which 
to integrate EHH. In such a case the demands of “compu¬ 
tational causality” and traditional physical/temporal causality 
differ, as has been discussed by Kinsler im (also see refer¬ 
ences therein). 
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dfv = apD - ypdiP - C(}|P 
: Vx [;U-'(B-M)] +J(0 



5,B = -Vx [e-'D-P] 
dfM = a^B - 7M'5rM - to^M 


FIG. 12: Causal diagram for the macroscopic Maxwell’s equations, 
with a unit cell being indicated by the dotted oval. Although fields E 
and H, have been added for clarity, with the arrows from B, M com¬ 
bining to indicate how the curl of H causes changes to D; and those 
from D, P combining to indicate how the curl of E causes changes to 
B. No parameters are added to the diagram to avoid clutter, but they 
are present in the equations. The locations where the parameters 
should be added, if so desired, can be inferred from the Lorentzian 
diagram in fig. and Maxwell’s equations. 



FIG. 13: Causal diagram for the TDDE. 

VI. OTHER EXAMPLES, WITH COMPLICATIONS 


Here Q is a source term, such as a driving term or some mod¬ 
ification to the wave equation. Next, 7] is a positive constant 
that imparts loss, and c is the high frequency speed of sound. 
The causal diagram for this is shown on fig. where the 
usual equation has been rearranged to conform to the explic¬ 
itly causal form with the highest order time derivative on the 
left hand side. 


B. An elastic rod wave equation 

Some mathematical models of physical processes, even 
rather established ones, are less easy to cast as a causal 
diagram. To make this point, I will consider an existing 
model representing acoustic waves traveling along an infinite, 
isotropic and elastic cylindrical rod of radius R. Following 
Murnaghan’s free energy model, Porubov has derived a wave 
equation governing propagation of the solitary waves along 
such a rod ca, and the propagation and dispersive proper¬ 
ties of this model in time propagated and space propagated 
pictures have been compared by Kinsler cni. There is no 
impediment in this model against the rod having “auxetic” 
parameters ifThl . e.g. where the Poisson’s ratio was negative 
ca. This “elastic rod equation” (ERE) describes the dis¬ 
placement g = g(x,t) with a second order wave equation of 
the form 

c^d^g - d}g + b^dfdlg - b2d^g + Xd^g^ = Q- (6.2) 

where g is the longitudinal displacement in the rod*. Here, 
when we move the highest order time derivatives to the left, 
we get 

(1 + bidl) g = c^dlg - b2d^g + xd^g^ - Q- (6.3) 

Most notably for my purposes here, this model has a prob¬ 
lematic feature: there are multiple terms in eqn. ( |6.2[ ) that 
have the highest (second) order time derivatives. These either 
have to be combined into a single term, to give us a single 
quantity that suffers the effects of the others, or the equations 
must be partitioned into two parts each of which is causal in 
appearance. However, here there is no straightforward way 
of achieving this here, because of the presence of the bid^ 
component; in simpler cases such as (e.g.) the E-model for 
the spilt ring resonator, it can be done (see sec. IV of ifTl). 


A. The time-dependent diffusion equation 

The time-dependent diffusion equation (TDDE) ifT^ is a 
second order wave equation with a loss term added; it appears 
in a variety of contexts in physics, including acoustic waves 
in plasmas or the interstitial gas filling a porous, statistically 
isotropic, perfectly rigid solid 03. It has a three-dimensional, 
inhomogeneous form for the velocity potential g = g{f,t) of 
na, and is typically written as 

V^g-c^^d^g-r]dtg = Q. ( 6 . 1 ) 


^ For the sake of completeness, the other parameters (following QSl) are: 


p 


E B , v(v-1)r2 v£R2 

Po 2po ’ ‘ 2 2po ’ 

3£: + / (/ - 2v)^ -f Am (I - 2v) (/ + v) + 6fiv^ 


Here jS is the nonlinear coefficient, E and v are Young’s modulus and Pois¬ 
son’s ratio respectively, /, m, n specify Mumaghan’s modulus, and po de¬ 
notes the density. Poisson’s ratio is typically rather small (i.e. |v| < 1), in 
which case b\ will be negative. In contrast bi can cover a wide range of 
values, especially if auxetic materials are considered, but usually v,E > 0, 
so that ^2 < 0- 
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Such terms also occur in other systems, such as the Van Wi- 
jngarden’s equation ifTOllTSl for waves in bubbly liquids: i.e., 
terms combining derivatives which are both second order time 
and second order in space (here, b\dfd^g). Other terms, such 
as the d,V^g terms appearing in both the Stokes’s m and 
Van Wijngarden’s equation do not cause difficulties because 
they are only first order in the time derivative, which is not the 
highest order present. 

The best strategy here would be to remrn to the original 
derivation of the ERE, and alter it: perhaps, for example, there 
might be an equation for how d^g responds to stimuli, so that 
eqn. could be split up in some physically motivated way. 
Eailing that, there are two alternate strategies we might use 

First, we might decide to use a spatial spectrum or “k- 
space” picture instead, replacing the space and time displace¬ 
ment g = g{x,t) with g = g{k,t). As part of the conversion, 
each instance of dxg then becomes a factor of ikg-, the result¬ 
ing 1 +bik^ factor on the left hand side can then be divided 
out. 


d}g{k,t) 


-ch^-b2k^ 

l+bik^ 


g{k,t) 


Q 

l+bik^' 


Xk^[g{k,t)*g{k,t)] 

l+bik^ 

(6.4) 


This expression can be diagrammed fairly easily, the only 
drawback being the more complicated calculations required 
to condition how g affects itself Also, that while this spa¬ 
tially non-local k-space representation can remain temporally 
causal, it does not respect finite signal speeds IfTTl . 

Second, we could reject strict temporal causality, and, as is 
often done in optics and acoustics nniini, replace it with 
what might be called a “spatial causality” instead. This is 
where we propagate our quantities of interest forward along 
some trajectory in space, rather than forward in time. Thus, 
instead of following how a spatial profile g{x) evolves as time 
t passes, we instead model how the time history g{t) evolves 
as we walk along the x axis. Since eqn. has one spa¬ 

tial dimension we can ignore the complications of transverse 
spatial behaviour, and treat x as if it were the time-like dimen¬ 
sion. Writing T = x/c and ^ = ct, we then proceed as if T 
really was the time axis. The relevant highest order derivative 
is the b 2 dxg one, but it gets converted to {b 2 /c^)d^g. With 
the rest of the conversion complete we find that the scaled 
time history g {^) follows 





c ^2 I ^ 




b2 

(6.5) 


Unlike the spatial spectrum approach given first, each of the 


cause terms generated here is relatively simple, however there 
are five of them rather than just three. 

In summary, the attempt to diagram this ERE model in a 
causal way has raised an important issue: not all model equa¬ 
tions have rigorously causal interpretations. This is perhaps 
not unsurprising, since many simplified descriptions of physi¬ 
cal phenomena are far removed from the starting point of their 
derivation, and the effect of simplifying approximations can 
be subtle and unexpected. 

VII. CONCLUSION 

One of the primary benefits of this technique is that in mak¬ 
ing the diagrams, we are forced to clarify what we mean 
by “changes”. It is not uncommon to hear, in conversations 
about EM, statements along the lines of “Maxwell’s equations 
tell us that changes in E cause changes in B”. Here we are 
forced to distinguish between changes in time (codified here 
as “effects”) and changes in space (codified here as potential 
“causes”). Generally, it is best to restrict the use of the word 
“change” to denote “effects”, i.e. changes in time. 

What we can now say, within a rigorously defined prescrip¬ 
tion Ul, “spatial variation in E causes changes in B”, or “spa¬ 
tial variation in H causes changes in D”, and even draw a di¬ 
agram to emphasize the point. Eurther, the diagram can leads 
us to an algorithm for computation in which calculation of ef¬ 
fects (changes in the state of the system) are driven by causes 
(previously calculated states of the system). 

It remains to be seen whether or not causal diagrams like 
these find a useful place in building or analysing physical 
models. Diagrams such as Eeynman ll20l or Wyld diagrams 
Il2ni22l have their fields of application in field theory or fluid 
mechanics, as a means of bookkeeping and calculation on 
perturbation expansions. Others, like Block diagrams lU or 
DAGS 0, and as already mentioned, are not designed with 
dynamical (time dependent) processes in mind. 

Here, however, the intent is for a more general type of dia¬ 
gram based simply on representing differential equations in a 
causal way, and as a side effect, also leading the way towards 
an algorithm for numerical simulation. As part of this, they 
also provide a way to reinforce a rigorously grounded notion 
of causality in even the simplest the physical models. Even 
F = ma, trivially rearranged to read dt\ = F/m has a causal 
diagram, one that matches the very first one presented here: 
the left hand panel of fig. [T] 
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